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SUMMARY 


A frequency-response method for determining the critical control- 
gearing and. hunting oscillations of airplanes wiili autcmatic pilots 
is pi'essnted. The method is grajhical and has several advantages 
over the standard nmerical procedure based on Rouih's discriminant. 
The chief advantage of the method is that direct use. can be made of 
the measured response characteristics of the autcaaatic pilot . This 
feature is especially useful in determining the existence', amplitude, 
and frequency of the hunting oscillaticns that may be present vhen 
the automatic pilot has nonlinear dynamic characteristics. 

Several examples are vorked out to Illustrate the application 
of the frequency-response method in determining the effect of 
automatic -pilot lag or lead on critical control gearing and in 
determining the amplitude md frequency of hunting. It is shown 
that the method may be applied to the case of a control geared to 
airplane motions about two axes. 


D^TRODUCTICN 


The increased rise of automatic control on aircraft (especially 
on pilotless aircraft) has focused attention on the proper design 
of the aircraft and associated control systems with a vlev toward 
obtaining satisfactory dynamic stability. The control system 
(autopilot) consists of the gyro, phase -shifting device, and 
servomotor. Factors in the control system that determine stability 
characteristics are control gearing, lag in the serycmotor, and 
lead in the phase-shifting device. The stability characteristics 
of the airplane depend on airplane ccnfiguration and mass distribution 
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The purpose of Idie present • report Is to give a method for 
analyzing the dynamic stability of an airplane with automatic 
control which mahes direct use of the observed dynamic charac- 
teristics of the automatic pilot (lierein ceilled autopilot) . The 
method separates the characteristics of the autopilot frcm those 
of the airplane. It, therefore, easily reveals the effects of 
modifications to the autopilot such as adding lead. The procedure 
is largely graphical. The method Is similar In certain respects 
to that of Kyquist (reference l) , which was devised for electronic 
circuits but which has also been applied by some workers to the 
design of servomotors. 

The frequency-response method was previously applied by 
Jones (reference 2) to calculate the hunting produced in airplanes 
with "flicker," that is, on-off control. If the • servomotor 
(herein called servo) has a linear lag characteristic' (lag indepKident 
of amplitude), xhen there will .usually be an upper limit to the 
control gearing above' ■vdxlch the airplane will- be unstable . The 
frequency-response method de.termlnes this crlticej. control, gearing 
and. the corresponding frequwcy for any type of lag or lead in the 
autopilot. It also indicates the changes to be made in the servo 
which will improve the' stability of the airplane. 

The utility of the frequency-response method is most appsirent 
in the case of an actual servo with nonlinear lag characteristics. 

This method, allows the use of measured characteristics of the servo 
that might .be inconvenient to represent by a mathematical formula. 


TEEMiNOLOG-Y AND SYMBOUS 


The word "lead" is tised In this report In two ways. It is the 
phase angle of control-surface deflection 6 ahead of airplane 
deflection, say angle of pitch 0, and it is also used to indicate 
a device that causes an increase in the phase angle of lead of 6 
ahead of '©. These phase-shifting devices are herein called 
first- derivative and second-derivative lead. 

The expression- "linear autopilot" is used to indicate a servo 
■that is acted upon b'y a force proportional to the input signal 
(airplane deflection) resisted by a force proportional to "the 
displacement and velocity of the ou'bput (cantrol-surface motion) . 

The function of -the automatic pilot is ' ■to apply a corrective 
control deflection in response to any deflection of ■the airplane. 

The "control gearing" is ■the ratio of the applied control deflection 
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to the airplane deflection for very slow deflections (static ocnditicn) 
There may exist a “critical control gearing' which in nearly all 
practical cases is an upper limit heyond which dynamic instability 
occurs. - . _ . 


R 

^cr 

kg 

k 

c 

I 

m 
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The following symbols are used: 

ratio of control deflection due to second derivative of 

airplane displacement to that due to airplane displacement 

ratio of control deflection due to first derivative- of airplane 
displacement to that due to airplane displacement 

amplitude .of airplane oscillation produced by unit amplitude 
oscillation of control surface 

amplitude of control- surface oscillation required to maintain 
unit amplitude of airplane oscillations (i/R) 

amplitude of control- surface oscillation produced by autopilot 
in response to airplane oscillation, divided by amplitude 
of airplane oscillation 

for unit control gearing 


critical value of control gearing 
spring stiffness factor of servo and control surface 
factor relating airplane deflection and force on servo 
control gearing, ratio of 8 to 0 in static condition 
viscous damping factor of servo and control 'surf aice 
viscous lag factor which depends on o and k^ (c/k3_) 


(^Ai) 


mass of movable part of servo and control surf^e 
time, seconds. 

period of oscillation, seconds 

angle of control- surface deflection in direction to reduce 
angle of pitch of aiz’plane, 'degrees 


8 



k NACA 3N No. 1229 

g phase angle of lead of- 8 ahead of 6 vhen oscillating 

^ airplane forces control surface to oscillate, degrees 

€ phase angle of lead of 6 ahead of B when oscillating • 

control surface forces airplane to oscillate, degi'ees 

9 angle of pitch of airplane, degrees 

8 signal fed Into servo 

^ angle of hank of airplane, degrees 

\|; ■ angle of yaw of airplane, degrees 

$ gyro displacement produced hy airplane deflectlcn 

05 angular frequency, radians per seccnd (2it x Frequency) 

05 angular frequency of servo and airplane when k » k^^,, 

radians per second . .. 

05j^ natural angular frequency of servo, radians per second 

T angle of tilt of gyro axis from X-axLs of airplane, degrees 

D differential operator (d/dt) 

Subscript; 

max maximum 

BASIC PEINCIPLES OP PREQUMCY-EESPCTfSE METHOD 


The response of the airplane to a sinusoidal control motion may 
be computed from the eqxiatlons of motion of the airplane . (See 
appendix A.) The airplane response is sinusoidal and of the same 
frequency as the control motion and has an amplitude and phase that 
depends on this frequency. If the ccntrol motion is given by 

■& » sin 05t 

then the response in pitch, for example, inay be expressed as 

© = E sin (cDt " €j,) 
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where R and both depend on m, and E Is the ratio of Idle 

amplitude of Q to the amplitude of 8 when the control surface is 
forcing the airplane to oscillate. Changing the fniDplitude of 8 
produces a proportionate change in the amplitude of 0 hut does not 
affect the phase angle fij,. The reciprocal of R gives the 

amplitude of 8 required to sustain unit ^plitude oscillation 
in "■ S and is denoted hy that la, if 

. 8 » sin cut ' . , 

then the response will he 

0 = sin (oit - €|p) 

A plot of and Cj, agdlnst- oi ‘is then made and Is referred 'to 
as th'e frequency response ‘of ’ the airplane. ' ’’ ! 

The frequency response pf the autopilot , is obtained 1?j5r measiirln.g 
or cckpu^ing the response of the’ ccntrol surface h/o a steady os’c illation 
of the airplane. The response may he calculated if the dynamic 
constants, of: the autopilot are taacwn, hut it is iisually ' easier to ,■ 
measure the response "by oscillating the autopilot in the laboratory. 
These calculations or measurements give the respcnse "to ' 


6 = sin cot 


and this response is ezpressed hy ■ . ..:.i 

8 = Kp^ sin '(cuf + 6p) 

where KL *ls the amplitude of .8 --for unit control geaning and €_ 

- ^ - .^1 - -. v - ; . - --- *: 

is the phase, pngle of lead of ■ 8 ahead of 0 when the airplane is 
forcing, the control surface to oscillate. 


If. at sc?ae- value of the -angular frequency co the condition 


exists^ then "the airplane comected to the autopilot will oscillate 
continuously at -this frequency at consist amplitude^ provided the 
control gearing is made equal to ’ at this v^ue of ojj 
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that is, the critical control geeirlng is given "by 


cr 


- 5 : 


at 0 ) where « £p. This ccaiditlon of neu-feral stahlllty is shown 
schematically in figure 1. Aa exemple of a solution for lc-„ is 


shown in figure 2. In this case o> 


cr 


6 .2 and k 


cr 


K_ 


at this 


value of oij therefore 


k _ 242a 

cr - 0*91 


0.36. 


In almost all cases control gearings greater than result 

in Instahility and ccxitrol gearings less than k^j, provide stability* 
Thus, kgj. is usually an upper limit to the permissible control 

gearing. More generally, the stable side of the boundary is determined 
by the relative slopes of the Cj, and c,, curves at their 

intersection. If is greater than the system is stable 

din do) 

below k^j, and vice versa. A proof, of this rxile is given in 
appendix B. . .. . . 

In the preparation of figure 2 and the examples that follow 
(figs. 3 to 8) a particular airplane configuration was assumed for 
the computation of and (See appendix A for methods.) 

This hypothetical pilotless airplane has four fine 90? apart, a 
horizontal pair for pull-ups and a vertical pair for turns. The 
fins constitute the only lifting surfaces of the airplane and are 
equipped with ailercns connected to a roll- gyro and servo that maintain 
the airplane in a position in which one pair of fins is always 
horizontal. The pitch and yaw control surfaces are tralling-edge flaps 
on the fins and are connected in pairs, one pair on the horizontal 
fins connected to a pitch gyro and servo and the other pair cn the 
vertical fins connected to a yaw gyro and servo. It is assumed that 
the pitch £aid yaw motions do not interact with each other ctr with 
any accidental rolling motion that may exist. Figures 2 to 6 show 
the elevator-pitch stability. Budder-yaw stability may be expected 
to be very nearly the same as elevator-pitch stability because of 
the conf igxiration of the fine. 
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Linear Autopilot 

The motion of 'a control surface actuated hy a simple linear 
autopilot is similar to that of a mass-spring dashpot and may he 
represented hy the oquatioh 

mD®6 + cD6 + Is^S = kgS . (l) 

•where each of ■the factors m, c, and is the sum of two parts, 

one part due to the ser'vo and one part due to ‘the control surface . 
For example, lC 2 _ is -the sum of the spring constant of -the servo and 

■the aerodynamic hinge-moment constant of the control surface . The 
■term represents a force proportional "to "the deviation in angle 

of pitch Q applied, -without lag. 

Equation (l) may he expressed as 



+ .£- D + JO. 

\ ^ ) 




where IsQ^/'kr^ is the con-trol gearing k. The natural frequency of 
the servo is given hy , _ . . 



and the damping depends on "^e quantity I » -2-, With these 

^1 

suhstitutions, equation (2) hecomes 




ke 


The term 7D represents .viscous damping and the term 
represen-ts inertia reaction. 



When the autopilot is oscillated in pitch, 'the control 
response 6 lags behind 9. The lag usually decreases as' -the lag 
factor I decreases and as the natxiral frequency increases. 
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It Is convenient to choose o:^ equal to 2/?, a value that gives 
critical damping. With this choice, the phase a n gle gp and the 

relative amplitude for unit-control gearing Kp of the response 

are presented In figure 3 four values of 2. Values of 
and €j. for the assumed airplane, are also given so that 
can he obtained for each value of I- 

If first-derivative lead is added to the system, the equation 
of motion relating 6 and 0 • becomes 


5 


(• 


m + 



= 3r(l + rf>)0 


Figure 4- shows E and for 1±ie same values of I and tuu 

Pi P "-n 

as figure 3; but vlth r *= i. Comparison s^ove that- the first- 

derivative lead produces an angle of lead that approaches 90° at 
high values of the angular frequency to. This amount of lead is 
enou^ to override the lag for moderate values of I at low 
frequencies but not at high frequencies. 


Values of obtained from figures 3 and 4 are shown in 

figure 5 plotted against the viscous lag factor Z. Ccneiderable 
increase in the stable range of control gearing is evident for 
small values of Z with the first-derivative lead. Note how -this 
Improvement in stablliiy diminishes at large values of Z» 


If both first- and seocaxcL-derivative lead are added, the 
response of the control is given by the egwtion 


5 


k 


1 rP 
1 + ZP + 



0 


It can be shown that If 


and 


r> I 


> f ■ 


(5) 
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lead will . exist at all' f xegutencies' . Iii this situation •^e ’airpla^eV;.'^ 
cannot .oscillate regardless of •the', value of k or^ ste'ted ' • ■ ■ 

dif-fer.ently', the critical control gearing is Infinite. Figure 6 
shows the resultant lag or lead for a critically damped "servo for 
•which oJq = 20 radians per second and for various ■values of r and a. 

The •top curve is for values of r ar.d a that satisfy expression (3) 
and, -therefore, result in lead over -the entire frequency range. 


Nonlinear Au-topllot 

Even if the servo is constructed -to give a response proportional 
to a steady disturbance, i-ts response -to an oscillating disturbance 
is not proportional in practice -to the amplitude of -the disturbance 
because of the nonlinear oynamic characteristics of -the servo. Thp. 
lag is then a function of bo-fch amplitude and frequency. By finding" 
experimentally -the frequency response of the autopilot for a nmber 
of amplitudes and for a given control gearing, it is possible -to 
determine the amplitude and frequency at which the airplane •will 
oscillate -when coupled -bo -the autopilot at -that particular value of 
control gearing. The condition for s-teady oscillatione (hunting) 
is that Ep = and = Cj. at same frequency. The frequency 

and amplitude of the s-beady oscillations of ■the airplane are -the 
values of tu and "which Ep '«? ■ and 

The stability of -these cons-bant-ampli-tude oscillations has a 
slightly -different significance from, that of ’ linear systms < ■ In 
linear systems 'the oscillations are, in gejaeral> ei^er damped or 
undampecL-regardless of -bhe amplitude. In--the"case of -the hunting 
•that may exist with a .nonlinear servo the , oscillations are s.ald -to 
be s-table if, after a dis-burbance from their s-beady value, they. -tend 
to return -to that steady value. The cri-berion for -the s-bability 
of the s-beady oscillations follows- from -tha-b for the line^ au-topilot 
and is " ' ... - 



• *1 . 

The solution ob-t^ined is, not exact, .because ■ -the . autopilot response 
to Bln-usqidal ai^^ans motion is no-t a pure .-sine wave but is dis-bortedj 
ho-vrever, the .eurfs tehee and aproroxima -be -amplitude and frequency of 
hunting oscillations of -the f'undamen-tal can be determined by the,;freguency‘- 
response -liethod-,- • Ther-valuds of- Ep and e-p ‘-used -when the response 
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is not sinusoidal are -tlie values of the equivalent sinusoidal response 
which supplies same energy per cycle and the seoae impulse per half 
cycle as the observed response . This condition leads to the following 
expressions for £p and 





max 


= tan"l S 
A 

■where 



The value of S in these Integrals is the measured control- surface 
deflection in response to an oscillating motion of ttie airplane 
represented by 0 » ®max O^rdlnarlly the oscillation in pitch 

is performed in the laboratory by mounting the autopilot on a table 
that can be made to oscillate . 

A typical application of the method for determining amplitude 
and freqiiency of hunting of an aircraft and autopilot with stops in 
the control system is shown, in figures 7 aod The airplane used 
in this example is slightly different from that considered previously 
and, therefore, the and curves ere slightly different froa 

those in flginres 2 to 6. Figure 7 contains the pheise and amplitude 
of the experimental response of the autopilot plotted against frequency 
for a series of amplitudes of the input and also the corresponding 
calculated curves for the aircraft. For each amplitude the values 
of Kp and Kj, are obtained at the frequency where fip •» The 

ratio Slid the value of (where $p » e^) are plotted 
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against amplitude of pitch in figure 8. The amplitiide and 

frequency of hunting occur vhere ~E. = i. In this example the 

amplitude of the hunting is about 2.5° and the frequency is 18.7 radians 
per second. . — . • . , . 

Nonlinearity In the autopilot usually causes an increase of -lag 
with amnlitude . This • increase may result in a reversal of the favor- 
able effect of first-derivative lead as the lead factor r is • 
indefinitely Increased, and suggests an optlrauni value for . r, a . - 

concluslco. that is supported by fli^t eaperience. This result may 
be traced to the counteracting effects of lead and amplitude intro- 
duced by the first-derivative lead. The first-derivative lead shifts 
the phase of the signal fed into the servo so that the signal leads 6 

by an amount tan“^ rca and increases the amplitude of the signal by 

• { 2 2 

the factor v 1 ♦ The change of lead with r may be stated 

symbolically as- follows: If the dependence of €p on 0 and to- ■ 

■without first-deri-vative lead is expressed by 


■then the addition of first-derivative lead changes the lead to 

€.p * f(8,co) +■ tan"^ rco -■ (k) 

where . ■ ■ 


6 = 0 \/l + r^cD^ 


Differentiating equation (4) 'wi-bh respect to r’ gives 


de 
— £ 
dr 


dc 
dD dr 


■ — ^0\/l""+~r^ii^- 


CO 


1 + 


d0 ■ 


CD 


'/l + .r^(D^ 


1 -t 




(5) 


If r is very small, the first term on ‘the right-hand side of equation (5) 

d£^ 

is negligible, so that --E. is positive. The addition, therefore, of 
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a small amount of first- derivative lead results in a change in €p 

in the direction of more lead. A similar analysis of the effect 

dK 

of r on K shows that — £ = 0 idien r = 0. For small values 
P dr 

of Vf thoreforo, the effect of adding first- derivative lead is the 
same as for the linear servo: an improvement In dynamic stability 

and a reduction in the amplitude of hunting is to be erpected. For 
large values of r^ however, the lag due to increase in amplitude 
(first terra) tends to co^mteract the direct lead (second tecnn) end 
may cause more violent hunting. Seme optimum value of r may 
therefore exist. A direct determination of-its value may be made 
from the autopilot for various values of r, Q, and to. 


Application of Method to Ccaatrols Geared 
to Motion about Two Axes 

Semetimes the aileron or rudder is geared to both the angle of 
bank and the angle of yaw. Usually this gearing is accomplishoi by 
tilting the gyro in such a way that the rotation of the gyro is 
affected by bank and yaw according to 1±io forraitla 

(5 ts ^ cos T + \}r sinT 

In order to apply the frequency-response method in this case it is 
necessary only to calculate the response in bank eind yaw of the 
airplane to a sinusoidal control motion and to combine them according 
to the preceding formula to obtain This formula yields the Ky 

and Cj, curves. The other pair of curves required - the Kp and 

curves - are obtained in exactly the same way as before, namely, by 
oscillating the gyro about its sensitive axis. The critical control 
geai’lng so deteimined will be the critical ratio between S and 


CCNCLUDIHG EESdABKB 


The frequency -resPCHise method of analysis is a useful graphical 
means of determining oscillation characteristics of an airplane 
equipped with an automatic pilot. If .the servo Is linear (ideal caee) 
the critical control gearing beyond which increasing oeclUatlons 
take place can bo readily determined. If the servo has nonlinear 
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characteristics (practical case), the- eiistence, amplitude, and 
frequency of steady hunting oscillations can he determined approxi- 
mately from the measured frequency response of the servo and the 
computed frequency response of the airplane. 

In general, the use of first-derivative lead (a phase -shifting 
device) has a fayorahle- effect on the dynamic stahility. Large 
amounts-. of first-derivative lead or phase shift may he destabilizing 
due to toe increase in toe signal amplitude produced hy this phase- 
shifting device . 


Langley Memorial Aeronautical Laboratory 

National Advisory Committee for Aeronautics 
Langley Field, Va. , December 9, 19^6 
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APPENDIX A 

CALOTLATICK OE RESPONSS OF M AIRPLANE TO SINUSOIDAL CCSSflROL MOTICH 


The response of the airplane in pitch to sinusoidal elevator 
motion will be; used as an example . The equations of motion, for a 
given control deflection, are— 



where 




Or, - ^ ^ 

^ pVS d7 


T = 

pVS 

% ” - - %) 



S =■ ‘x 

'a = 
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in 

__3L 
me 


and 

T 

m 

c 

S 


C « V -3 
“Si ^ 


Sc 

c = 

“a ^ 


trim angle of clinib, deg 
thrusty It) 

pitching moment of Inertia, eltig-ft^ 
mass of airplane, It 
wing chord, ft 
wing area, sq ft 


u « AV/T 

V velocity, ft/sec 

AV change in velocity, ft/eec 

lift coefficient 
a angle of attack, deg 

pitching-moment coefficient 
Cjj ‘ drag coefficient 

p air density, slugs/cu ft 
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Solving eqaiations (Al) for 

gives 


i. "by the method o£ determlTiantB 
S 


s 



+ tB 


0 


z 

u 

Za + T® 

-^5 



“a “Da 

“s 


+ TD 

a 



z 

u 

z„ + TD 
a 

Zq - TD 


C 

m 

u 

C + C D 
“Da 

C B - p.]j. 2j 

“ne 2v "^y ' 


Expanding these determinants hy the ■usiial methods gives an expression 

for £ consisting of the ratio of tvo polynomials in V, which may 
6 

he written as 


Q ' aj^I^ -*• h^D + c-^^ 

6 k •:} p 

a^D + hgD-^ + CgD^ + 6qD + ©2 

The phase and amplitude of the response of 0 to the motion 6 = sin cot 
is obtained hy substituting ixo for D In the above expression for 

This substitution gives a complex number, say, A + IB. The angle 

of lead of 0 ahead of ^ is tan'’^ S and the amplitude C3f 0 

is \/a^ + that is 

«r ' I 

R *= '\/A^ + B^ 

0 => E sin /cot - 

^ 

^ R 


and 


0>JQ> 
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APPMDIX B 

DETEHMiNATICN OF STABLE SIDE OF STABIUTS: BOUNDARY 
OBTAINED BY FBEQUENCY-KESRCKSE METEOD 


It is shown in the hody^ of this report that a critical value 
of the control gearing exists if at some frequency oj the veilues 
of Sj, and Sp are equal and that the value of the critical 



he shown timt the stable side of the stability boundary may be 
determined by the relative slopes of the £ 3 , and £p curves at the 
point where « Sp. - • - 

As shown in reference 2, the resTpanse of an airplane to any 
disturbance may be computed fraa the response function.' The 
response to a sinusoidal disturbance is obtained by substituting lui 
for the variable D in the response function. The response functions 
for the airplane and autopilot may- be cambined to fom the stability 
equation of the airplane plus autopilot as follows: 


■ Let 6 = f]_(D)5 be the airplane response function and 
S * fc f 2 (D )0 be the autopilot response function where k is the 
control gearing. Combining the two gives ' 

5 « k fgCD) fx(r >)6 


or 


1 - k fgCD) fi(D) = 0 

By substituting f (D) = f^Cl^) f2(®) this equation may be written as 

1 - k f (p) « 0 (Bl) 

This is the "stability equation," the roots of which for the variable D 
determine the damping and frequency of the motion. The frequency 
responses are obtained from the response functions by use of the 
following: 

f, (ico) 55 Re a J.- e ^ 

^ Kp . 

fgClcj) «s Kp^e ^ 
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Then a d^mamlc-stahlllty ‘boundary exlata if 1 - k f (la>) = 0 
for scane value of eay a»Qy* follows from the definition 

of f^_ and that this dynaialc-Btablllty boundary occurs when 



and 


Jr 


‘cr 


at some value of o) * Under these conditions, equation (Bl) has 

a pair of roots D » t IcCcy when k has the value given. 


In order to find the stable side of the boundary, it is necessary 
to find the sign of the real part of © at the boundary. From 

dk 

equation (Bl) 


dk 




kf *(D) 


This equation must be evaluated Tdien D <= At this veilue of ca, 

or ^ 

f(l£DQy) » :r~ • In genereil, 



Then 


f ' (p) « df(dfl>l ^ df(to) 
d(io>) d£0 


therefore 


^(l(of 

d(io>) 


<D=0), 


cr 
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end 



If 



is negative, the real part of 


o 

4 ^ is positive 

dk 


As k increases above the critical value, therefore, the roots of 
the stability equation indicate instability; that is, if 


4p - ^J) 


do} 


< 0 


_| 0 )=< 0 , 


'or 


then the stable region is located -vdiere k< k^j, and vice versa. 
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K 
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cr 


K 
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Figure llluetration of frequenojr-respoaee aetbod of deteraiiLlng oritloal 

control gearing. o>or * ^or ** (Dotted line Indicates where 

■ «j.. ) 
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Fig. 3 
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